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ABSTRACT 

Nova Cassiopeiae 1993 (V705 Cas) was an archetypical dust-forming nova. It displayed 
a deep minimum in the visual light curve, and spectroscopic evidence for carbon, 
hydrocarbon and silicate dust. We report the results of fitting the infrared spectral 
energy distribution with the dusty code, which we use to determine the properties 
and geometry of the emitting dust. The emission is well described as originating in 
a thin shell whose dust has a carbonisilicate ratio of 2:1 by number (~ 1.26:1 by 
mass) and a relatively flat size distribution. The 9.7 /zm and 18 /xm silicate features 
are consistent with freshly-condensed dust and, while the lower limit to the grain size 
distribution is not well constrained, the largest grains have dimensions ~ 0.06 /im; 
unless the grains in V705 Cas were anomalously small, the sizes of grains produced in 
nova eruptions may previously have been overestimated in novae with optically thick 
dust shells. Laboratory work by Grishko & Duley may provide clues to the apparently 
unique nature of nova UIR features. 

Key words: circumstellar matter - stars: individual: V705 Cas - novae: cataclysmic 
variables - infrared: stars 



1 INTRODUCTION 

V705 Cas (1993) was a typical dust-forming nova which 
showed evidence for carbon, silicate and hydrocarbon dust 
(Evans et al. 1997, Mason et al. 1998). Early lUE obser- 
vations (Shore et al. 1994) suggested that the grains grew 
to ~ 0.2 yLxn shortly after dust condensation. Mason et al. 
(1998) modelled broadband infrared (IR) data using a com- 
bination of dust types, and concluded that the circumstellar 
dust shell consisted of carbon and silicate grains in the ratio 
~ 8.5:1 by mass. 

Evans et al. (1997; hereafter Paper II) presented IR 
spectroscopy of the dust in V705 Cas, with some data cov- 
ering the wavelength range 2-24 pm. In addition to nebular 
and coronal features, emission features normally associated 
with hydrocarbons - the so-called 'Unidentified Infrared' 
(UIR) bands - were present in the spectra. The UIR fea- 
tures in V705 Cas were at wavelengths 3.28, 3.4, 8.1, 8.7 
and 11.4 ^m; the '8.1' and '11.4' feature seem to correspond 
respectively to the UIR features normally seen at 7.7 pim and 
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11.25 /im (but see §4.6.2 below), while the '3.4' feature was 
much stronger relative to the '3.28' than is generally the case 
in other astrophysical environments. The properties of the 
UIR bands in V705 Cas and in other novae indicate a set of 
UIR features peculiar to novae (Geballe 1997), due either to 
the nova-specific environment or excitation conditions. 

The 10 ^m spectrum of V7Q5 Cas also displayed a 
prominent silicate feature at ~ 10 /im and possibly a broad, 
weak feature at 18 /xm (Paper II) . The 10 /^m feature in 
V705 Cas was narrow, and peaked at 9.7 /xm, in contrast to 
the broader silicate features in novae Aql 1982 and Her 1991, 
which were also shifted to longer wavelength (see Smith et 
al. 1995). This implies that the silicate in novae Aql and Her 
displayed a degree of crystallinity, whereas that in V705 Cas 
was amorphous. 

In Paper II we fitted a simple function of the form 
B{iy,T), where B is the Planck function at temperature 
T and /3 (~ 1) is a constant (the so-called /3- index for the 
dust), to IR spectra in the range 2-24 /^m, and concluded 
that the grains eventually grew to ~ 0.7 /im. However, given 
the high quality of the data and the availability of the dusty 
code - which solves the radiative transfer in a dusty envi- 
ronment for a variety of conditions (Ivezic & Elitzur 1997, 
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Figure 1. (a) Fit of dusty model to 1994 August data for nova V705 Cas; thick curve is dusty model. Note the strong UIR features in 
the 10 /nm window; the apparent 'noise' in the ~ 2;um window is caused by the nebular lines, (b) As (a) but for 1994 October/November 
data. 



Ivezic, Nenkova & Elitzur 1999) - we are now in a position 
to make a more sophisticated attempt at fitting the data 
and to re-evaluate the nature of the dust. 

Here we describe the apphcation of dusty to model the 
dust around V705 Cas. 



2 THE DATA 

The data analysed here axe described in detail in Paper II; 
for consistency with earlier work we take reddening E{B — 
V) = 0.5 (Hauschildt et al. 1994) and 1993 December 14 as 

the origin of time. Wc revisit the distance below. 

The fitting procedure works best if there are data on 
either side of the peak of the dust emission. For V705 Cas, 
two datascts satisfy this criterion, namely those for 1994 Au- 
gust (to which we shall refer as Epoch 1), and 1994 Octo- 
ber/November (Epoch 2); at these two epochs we have quasi- 
simultaneous data over the wavelength range 2-4 fim {K, L), 
7.5-13 /xm (N), 16-24 pm (Q), at resolution ~ 300 - 1000 
{K,L bands) and ~ 60 {N,Q bands). The Epoch 1 data 
were obtained within 4 days of each other (on days 251 and 
255 of the outburst, from Paper II; we take t = 253 days), 
while the Epoch 2 data were obtained within 40 days of each 
other (days 300 and 341; we take t = 320 days). The data 
are shown in Fig. 1. 

In each case wc assume that the nova and its environ- 
ment did not change substantially between the times the 
data were obtained. Obviously this assumption is less se- 
cure for the October/November data but our attempts in 
Paper II to fit the data with the simple function as discussed 
above suggest that this assumption is not unreasonable. 



3 THE FITTING 

3.1 The state of V705 Cas in 1994 

3.1.1 Preamble 

We use observations of V705 Cas to estimate the likely prop- 
erties of the nova and its environment at the time of our 
observations. For the purposes of the dusty fitting we shall 



simply assume that the stellar remnant radiates like a black- 
body. While more sophisticated models of the nova pho- 
tosphere are now routinely available (see Hauschildt et al. 
2002 and references therein), we have no a priori informa- 
tion about the state of ionization of the ejecta; however we 
expect that they are transparent to photons in the Lyman 
continuum at the times of interest. 

The effective temperature T, of a nova radiating at con- 
stant bolometric luminosity rises as the pseudophotosphere 
collapses back onto the white dwarf (Bath & Shaviv 1976). 
Consequently the emission shifts to shorter wavelengths and 
the visual flux (in mag) declines according to 



|Amvis| ^ 7.5 logio (^) 



(1) 



where Arrivis is the decline in mag from maximum light and 
To ~ 15,280 K (Bath & Shaviv 1976). However, observa- 
tional evidence is that novae at maximum are on the whole 
rather cooler than 15,280 K, a value closer to ~ 8, 000 K 
being more typical (Beck et al., 1995; Warner 1995); indeed 
the lower value is very close to the effective temperature de- 
termined for V705 Cas around maximum by Hauschildt et 
al. (1994). We take To = 8, 000 K here and use Equation (1) 
to estimate T*. 

We have taken the light curve in Paper II (reproduced 
in Fig. 2a) and dotorminod rhvis by a least squares fit over 
the period t = days to t = 60 days, omitting 'glitches' such 
as the small pre-dust 'dip' noted in Paper II. We find 



mvis = 6.32(±0.04) + 0.0364(±0.0012) t , 



(2) 



where t is in days; we therefore have that rhvis = 0.0364 ± 
0.0012 mag day"^ {ta = 82.4±2.7 days, t2 = 55.0±3.6 days, 
tn being the time for the visual light curve to decline by 
n mag), which differs somewhat from that assumed in Pa- 
per II (rhvis = 0.044 mag day^^, ts ~ 68 days) and is more 
in line with that given by Mason et al. (1998). We take 
^3 = 82.4 days here. 

Using the usual relationship between absolute magni- 
tude and rate of decline (Delia Valle & Livio 1995), the 
absolute visual magnitude of V705 Cas is My = —7.05 at 
maximum light; using a bolometric correction correspond- 
ing to an 8,000 K supergiant the corresponding bolometric 
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Days from maximum t (days) 

Figure 2. (a) Visual light curve of V705 Cas (filled squares), with Epochs 1 and 2 indicated; open circles arc estimated location 
of pseudophotospheric continuum, based on DUSTY fits and after applying interstellar extinction corresponding to E{B — V) = 0.5. 
(b) Time-dependence of the optical depth in the dust shell during light curve recovery, estimated as described in text. This is effectively 
the visual light curve of V705 Cas but rectified by Equation (2), which is represented by the horizontal line. The curve is an estimate of 
the dust optical depth after the deep minimum; open circles are values of r from dusty fits. See text for details. 



luminosity is 5.6 x 10* Lq. The value of Mv, together with 
the value of rrivisit ~ 0) from Equation (2) and the assumed 
reddening, gives a distance of D = 2.3 ± 0.4 kpc (cf. 3 kpc 
assumed in Paper II), in which the uncertainty in the red- 
dening is the major contributor to the error. 

3.1.2 Effective temperature of the remnant 

We estimate the effective temperature of the stellar remnant 
using Equation (1). For Epoch 1, we find T. ~ 135,000 K, 
while for Epoch 2, T, ~ 285, 000 K; there is a - 10% uncer- 
tainty in T* that follows from the uncertainty in the coef- 
ficient of t in Equation (2). Although there are in principle 
several ways of corroborating these values, complementary 
data at other wavelengths are generally lacking. 

In Paper II we noted the presence of [Sivi] 
(^Pi/i—^Ps/i) at A = 1.963 fj,m, which was fairly strong at 
Epoch 1 and also present at Epoch 2 ([Si vii] at A = 2.48 /j,m 
was also present, but very weak, at Epoch 2). The ionization 
potential of Si''+ is ~ 167 eV. If we assume that the gas is 
photoionized rather than coUisionally ionized (see Benjamin 
& Dinerstein (1990) and Evans et al. (2003) for a discussion 
of this), the presence of the [Sivi] feature implies a suffi- 
cient number of photons with energy > 167 eV. For strong 
[Si Vi] emission we suppose for the purpose of estimate that 
half the emitted photons should be capable of ionizing Si*"*" , 
implying an effective temperature T, ~ 1.93 x 10® K. At 
T, ~ 285, 000 K, only ~ 2% of the emitted photons are ca- 
pable of ionizing Si''+. The presence of the [Sivi] coronal 
line therefore seems to imply a T* that is very much greater 
than 285,000 K; alternatively, a different ionization mecha- 
nism may bo operating. 

On the other hand, Orio, Covington & Ogclman (2001) 
have looked at ROSAT data for over 100 classical and re- 
current novae. Although V705 Cas is not included in their 
discussion, they conclude that the 'super-soft x-ray' phase of 
the nova outburst is relatively short-lived, and is observed 
for only ^ 20% of novae: hot post-nova white dwarfs are 
relatively rare. Orio et al. hence estimate a general Icr up- 
per limit T* < 300, 000 K; the x-ray evidence (such as it is) 



seems to point to a value of T, consistent with that implied 
by Equation (1), and somewhat lower that that implied by 
the coronal lines. 

In our preliminary investigations (Evans et al. 2002) we 
approximated the star by a blackbody with T« = 70, 000 K, 
and noted that the model IR spectrum is not sensitive to 
the temperature of the central object. However the predicted 
ffux at A = 0.55 /im must not, after allowing for interstellar 
extinction, exceed the value implied by the visual light curve. 
This requires that the calculated by dusty must be 

^11, which provides a further constraint on T«. 

3.1.3 Dust shell thickness and optical depth 

CO was present in the IR spectrum of V705 Cas before max- 
imum (Evans et al. 1996, hereafter Paper I). In Paper I we 
assumed that the relative thickness of the CO-bearing shell, 
AR/R (where R is the inner shell radius), was ~ 0.1 in 
V705 Cas. Furthermore consideration of the dust-forming 
potential of novae, together with the rapid recovery of the 
visual light curve after deep minimum, supports the view 
that the dust too is confined to a thin shell (Rawlings & 
Evans 2005). The value of AR/R for the dust shell is pro- 
vided by our fitting procedure and, as the chemistry that 
leads to grain nucleation and grain formation is expected 
to be confined to the CO-bearing region, the dusty fitting 
may provide a posteriori justification for the assumption 
AR/R 0.1. 

Dust precipitation began on day 62, after which the vi- 
sual light curve went into deep decline. At this time, the 
difference between the observed visual magnitude and that 
predicted by Equation (1) is a measure of the visual optical 
depth r of the dust shell; the time-dependence of r is plotted 
in Fig. 2b. This is the visual light curve of Paper II, recti- 
fied by Equation (2), so that the visual decline described by 
Equation (2) is represented by the horizontal line. As the 
dust shell disperses (or breaks up) and the light curve re- 
covers, the pseudophotospheric emission shifts to the ultra^ 
violet (UV); at this time most of the visual light is provided 
by strong emission lines and Equation (1) no longer applies; 



4 A. Evans et al. 



Table 1. Nova and dust shell parameters for V705 Cas for the two epochs. See text for discussion. 



Property 


Epoch 1 


Epoch 2 


t (days) 


253 


320 


F(AC) 


64+° i° 


70+°- IS 


F(Silicate) 


o^+o.io 

'J-3O_0.04 




Td (K) 


604_+^6 


612112 


amin (Mm) 


0.003 ± 0.013* 


0.005 ± 0.025* 


Clmax (/ini) 


0.061 ±0.008* 


0.055 ±0.018* 


q 


2.3 ± 0.5* 


2.3 ±0.5* 


i?in (cm) 


7.181° ji X 10i« 


7.36l?;^s X 1015 


^in (mas) 


417.31^^7 


427.81^,^.6 


1 + AR/R 


1.11 ± 0.04* 


1.24 ±0.19* 


T 


0.0851°:°?^ 


19+0-20 


h.3 (10-1° Wm-2) 


3.26 


3.26 


/dust (10-12 Wm-2) 


13.1 


4.34 


/dust/ /2.3 


0.040 


0.013 


T* (K) 


135,000 


285,000 


Estimated Vphot 


12.8 


14.8 



for this reason the behaviour of the optical depth 'data' in 
Fig. 2b for t ^ 165 days is not meaningful. 

We see from Fig. 2 that maximum dust optical depth 
in the visual, Tmax — 5.5, was attained at t ~ 90 days, and 
we estimate that the optical depth for Epoch 1 was ~ 0.1, 
for Epoch 2 ~ 0.05. The precise values are not important 
at this stage; however our fitting procedure determines the 
optical depth in the dust shell and our analysis of the data 
using DUSTY should produce values of this order. 

S.l.J^ Dust composition 

The presence of the UIR hydrocarbon features is a strong 
indication that a form of carbonaceous dust is present. We 
therefore assume that the carbon component is amorphous 
carbon (AC), for which the optical constants are taken from 
Harmcr (1988). As indicated above, there is also a silicate 
component, and we take optical constants for warm silicate 
from Ossenkopf, Henning & Mathis (1992). 

We use our fitting procedure to determine the relative 
amounts of amorphous carbon and silicate dust in the shell. 

3.2 Fitting the dusty models 

For a given central source effective temperature, dust com- 
position and dust shell properties, dusty calculates the rar 



diative transfer through a spherically symmetric dust shell 
and determines the observed spectral energy distribution 
(SED). dusty exploits the 'scale-independence' of the ra- 
diative transfer problem, so that the observed SED is efi:ec- 
tively determined only by optical depth in the dust shell; 
this moans that 'absolute' values (e.g. luminosity and dust 
shell dimensions) are not uniquely determined by the trans- 
fer problem and must be inferred by other means. We use a 
downhill simplex routine (see Press et aJ. 1992) to explore 
the multi-dimensional parameter space and to determine the 
best fit between the data and the DUSTY output (see Tyne 
et al. 2002 for details). 

We assume an AC/warm silicate mix with a dust 
density distribution, appropriate for a steady wind. We op- 
timize the fit by varying the carbon:sihcatc ratio, the dust 
temperature at the inner edge of the dust shell, the grain size 
distribution, and the geometric and (visual) optical thick- 
ness of the dust shell. Spectral features, namely the UIR 
features, and nebular and coronal lines, were removed be- 
fore fitting the dusty output to the data, and restored once 
the fit had been optimized. The values of t, /g.s (the bolo- 
metric flux of the stellar remnant normalized to 2.3 kpc) and 
T, for each epoch were fixed. 

Further, as the relative placement of the spectra iu the 
K, L, N and Q bands is uncertain by ~ 10 — 20% (a conse- 
quence of the uncertainty in the flux calibration), the relar 
tive placements of the individual spectral components were 
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adjusted in an iterative way to optimize the fit. For each 
KLNQ dataset, dusty was called typically about 50-60 
times per optimization; including the iteration of the in- 
dividual KLNQ bands, dusty ran some 300-400 times for 
each epoch to optimize the fit. 



4 RESULTS &: DISCUSSION 

The best fits to the data for the two Epochs arc shown in 
Fig. la,b; these are equivalent to Figure 3a,b of Paper II ex- 
cept that in the latter, there was no attempt to optimize the 
band-to-band fit (sec §3.2), and the peak of the '8.1' band 
was omitted. The best fit parameters are listed in Table 1, in 
which F{AC) and -F(Silicate) are respectively the fractions 
(by number) of amorphous carbon and silicate grains in the 
dust shell; by definition, F{KG) + F(Silicate) = 1. Td is the 
dust temperature at the irmor edge of the dust shell, amin 
and amax are the minimum and mgiximum grain size in the 
grain size distribution 

n(a) da OC a da (ttmin ^ o ^ ttmax) ; 

here a is grain radius, n{a) da is the number of grains with 
radius £ [a,a + da], and 9 is a constant. 

In Table 1, /dust is the power output from the dust 
shell (at 2.3 kpc), and r is optical depth in the V band. 
For distance 2.3 kpc, the bolometric luminosity of the nova 
is 5.6 X 10* Lq at both Epochs (see Table 1 and §3.1.1), a 
consequence of our assumption of constant bolometric lumi- 
nosity (see Equation (1)). 

The uncertainties in a given parameter were obtained 
by fixing all other parameters and varying the given param- 
eter so as to increase the reduced by unity on either side 
of the minimum; where this procedure failed to converge in 
a reasonable number of iterations and only one error bound 
was obtained, this error is used and is indicated by an as- 
terisk in Table 1. The values of iiin follow from the input 
radiation field and fitted Td (see Ivezic et al. 1999); the un- 
certainties in iiin are derived from the corresponding errors 
in 6'in. 

Also given, in Table 1 and on Fig. 2, is the estimated 
contribution Vphot of the pseudophotosphere to the observed 
visual magnitude, based on the pseudophotospheric temper- 
ature, r and interstellar extinction; the estimated values of 
VpYiot clearly satisfy the requirement that they lie below the 
visual light curve (see §3.1.2). 

4.1 The dust temperature 

The dust temperature at the inner radius of the dust shell 
is determined to be Td = 604 K at Epoch 1, and 612 K at 
Epoch 2, i.e. essentially constant between the two Epochs 
within the uncertainties. These values are comparable to 
those determined by Mason et al. (1998, their Figure 1), on 
the basis of fitting Planck functions to broadband data for 
t ^ 275 days. However if the central star delivers constant 
power input to the inner radius of the dust shell, which 
consists of a single grain type that is unchanging (so that the 
grain properties are also unchanging) , we must have that 

. 2/(/3+4) 

Ta = n (f) , (3) 



for a dust shell moving away from the site of the explo- 
sion at uniform speed; here Tc is the grain condensation 
temperature and tc 62 days; sec Fig. 2a,b) is the con- 
densation time. From Equation (3), wc expect the value of 
Td t^'''-''^'*'' to be ~ constant, whereas for V705 Cas it rises 
from ~ 5524l5|g K day° '' (Epoch 1) to - 6149l^^i k day° '' 
(Epoch 2), assuming l3 = 1 (the value assumed is not crit- 
ical). This behaviour, in which the dust temperature does 
not decline in accord with Equation (3), is typical of dusty 
novae displaying 'isothermal' behaviour. 



4.2 The optical and geometrical depth 

The geometrical thickness of the dust shell, AR/R, is ~ 
0.11 ± 0.04 at Epoch 1, and ~ 0.24 ± 0.19 at Epoch 2 (see 
Table 1). In Paper I wc assumed that the relative thickness 
of the CO-bearing shell, AR/R (where R is the inner shell 
radius), was ~ 0.1 in V705 Cas. While this choice was some- 
what arbitrary it was governed by the fact that, in a nova 
wind that declines with time, there will be a thin dense shell 
at the outer edge of the cjecta (cf. Kwok 1983). Our DUSTY 
fitting indicates that the dust did indeed seem to have been 
confined to a thin {AR/R ~ 0.1) shell at the outer edge of 
the cjecta. Indeed, the similarity of the dust temperatures 
at the inner shell determined here (see §4.1), and those de- 
termined by Mason et al. (1998) on the assumption of an 
isothermal dust shell, is a consequence of the fact that the 
dust shell is geometrically thin so that the isothermal ap- 
proximation is a good one. 

Very roughly, for an optically thin dust shell in which 
grains scatter isotropically, the ratio of the dust flux to the 
bolometric stellar flux is ~ the optical depth through the 
shell, suitably averaged over wavelength. However the car- 
bon and silicate grains in V705 Cas do not scatter isotrop- 
icically and so this relationship will only be approximately 
correct; nevertheless it is gratifying to see in Table 1 that 

/dust 7/2.3 ~ T. 

4.3 Dimensions of the dust shell 

The inner radius of the dust shell, Rin, is 7.18 x lO'^^ cm 
(7.36 X 10^^ cm) at Epoch 1 (Epoch 2). At 2.3 kpc, the 
corresponding angular diameter 9 of the inner edge of the 
dust shell is 417 mas (427 mas) at Epoch 1 (Epoch 2). 

It is of interest to compare these values with the work 
of Diaz, Costa & Jatenco-Pereira (2001), who resolved the 
remnant of V705 Cas at 2.1 /xm at t = 2505 days (relative to 
the time origin assumed here). They determined an angular 
diameter of 840 mas, and concluded that emission by the 
resolved remnant they observed could not be free-free, and 
that the SED of the remnant could be represented by a 
blackbody at 1100 K. This value led Diaz et al. to conclude 
emission by dust; however our results, as well as those of 
Mason et al. (1998), show that the dust temperature was as 
low as ~ 600 K at t ~ 250 days, and it is difficult to see 
how the dust temperature could be elevated to 1100 K by 
t = 2505 days. 

ISO observations of V705 Cas (Salarna et al. 1999), at 
t = 950, 1265, and 1455 days, failed to detect any dust, with 
a limit of 0.3 Jy (~ 1.4 x 10"" W m"^ /um"^) at ~ 2.5 pm, 
the lower limit of the ISO SWS. The flux in the resolved 
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shell at 2.2^m was 4.1 x 10"^^ W m~^/xm~^ (Diaz et al. 
2001) so on this basis we can not rule out the possibility 
that these authors detected the dust shell. 

However if we and Diaz et al. have observed the same 
material, we would expect the angular diameter at Epoch 1 
(Epoch 2) to be 85 mas (107 mas) for uniform expansion, sig- 
nificantly smaller than the ~ 420 mas implied by our dusty 
fitting. There are two likely reasons for this discrepancy: 
First, if the emission seen by Diaz et al. is in the form of emis- 
sion linos (e.g. Pa-,5 A1.28 fan, Br-7 A2.16 ^m), the J and K 
band fluxes would be elevated relative to i/, £is observed by 
Diaz et al. For example, in data obtained for our programme 
on 1996 August 22 {t = 982 days), the P0 and Br7 fluxes arc 
1.7 X 10""' Wm"^ and 3.2 x 10"^^ Wm"^ respectively. Us- 
ing the filter bandwidtlis for the system used by Diaz et al., 
together with their J and K-hand fluxes for the extended 
sheU, their in-band fluxes are ~ 2.9 x 10~^^ Wm~^ (J) and 
~ 1.4 X 10~^^ Wm~^ (K) respectively. Given the difference 
in time between our 1996 August observation and that of 
Diaz et al., and the likely decline in the line fluxes between 
the observations, it seems plausible that they detected ex- 
tended emission from recombination lines. 

Second, with the Td values from Table 1 and bolometric 
luminosity L* = 5.6 x 10''Lq, blackbody grains would lie 
at distance ~ 7.5 x 10^^ cm from the stellar remnant, an 
order of magnitude smaller than the -Rin values in Table 1. 
On the other hand, spherical graphitic grains of radius Omax 
and having the Planck mean absorption efficiency given by 
Gilman (1974) would be at ~ 5 x lO^'^ cm. The deduced 
Rin depends on both the nature of the grain material and 
(via the T«-dependence of its Planck mean absorption ef- 
ficiency) on T*; the situation is further complicated in the 
case of V705 Gas by the grain mix and by the relatively 
weak constraint on T, (see § 3.1.2). This underlines the fact 
(see above) that the dust shell dimensions are not uniquely 
determined by the radiative transfer problem. 



4.4 The grain size distribution 

The grain size distribution n(a) da oc a~'' da, with q ~ 2.3, 
with no significant change in q between the two Epochs. In 
the environments of evolved stars, it is usual for the value of 
g to be ^ 3, a value typical of grain-grain shattering (Hellyer 
1970; Jones et al. 1995). The size distribution wo have de- 
duced for V705 Gas is significantly flatter than this. Such 
a situation may arise if grain growth/coagulation, destruc- 
tion and shattering in the wind of V705 Gas competed in 
such a way that the larger grains survived at the expense of 
smaller, possibly as a result of the former sweeping up the 
latter. Given that the charging of grains in nova winds has 
a significant impact on the rate at which grains grow, coag- 
ulate and ablate (Shore & Gehrz 2004; Rawlings & Evans 
2005), this is not unexpected; however we do not pursue this 
in detail here. 

Evans & Rawlings (1994) have concluded that free- 
flying polycyclic aromatic hydrocarbon (PAH) molecules 
would not survive the harsh nova radiation field. The sizes 
of the smallest grains in the distribution are not well- 
constrained, so we can not confidently determine whether 
the UIR features arise in free-flying PAHs or from the vi- 
bration of G— H bonds on the surface of hydrogenated amor- 
phous carbon (HAG) grains. 



The maximum grain size, on the other hand, is rather 
better constrained. We find that the largest grain size is 
~ 0.06 jum at both Epochs. This is significantly smaller than 
deduced in Paper II (^0.54/im at Epoch 1, >0.57/im at 
Epoch 2, and asymptotic dimensions ^ 0.7 /im) on the basis 
of \~'^ B(T, X) flts to the data; such discrepancies are not 
unexpected given the simple-minded treatment in Paper II. 
However it is also significantly smaller than the a^0.2/im 
(at t = 65.5 days) deduced - within days of dust precipi- 
tation - by Shore et al. (1994), on the basis of fiat UV ex- 
tinction shortward of 0.28 /xm, and is substantially smaller 
than the grain size generally deduced in nova dust shells, 
~ 0.5 — 1 /im (e.g. Gchrz et al. 1998). However we have per- 
formed several tests with large initial values of amax and 
flmax — 0.06 /im results irrespective of the starting value in 
the fitting routine, and we consider that this conclusion is 
robust. 

Wc believe that the discrepancy between our Omax and 
that deduced by Shore ct al. (1994) may have arisen, in part, 
for the following reason. Shore ct al.'s conclusion was based 
on the fact that the UV extinction was neutral shortward of 
A ~ 28OOA; they assumed that this implies 27ra/A > 4, so 
that a > 0.18 /im. While this is valid for primarily scatter- 
ing grains in the interstellar medium (see Spitzer 1978), it 
is less so for strongly absorbing (e.g. carbonaceous) grains 
in circumstellar environments, for which the condition for 
neutral extinction is closer to 2iTa/\ > 1. Indeed, the IR 
spectroscopic evidence is that the AG grains were the first 
to condense in the wind of V705 Gas (see §4.5), so the UV 
extinction seen by Shore et al. would have been due to AG 
rather than silicates. For V705 Gas therefore, the condition 
2na/\ > 1 is likely the more appropriate, so that neutral 
extinction implies a ^ 0.045 /tm. This would be more in line 
with the Omax values deduced here and would imply that, 
after an initial (~ 10 day) phase of very rapid increase in 
grain size, grain growth levelled off, possibly as the conden- 
sate was depleted. 

More generally, grain size in nova winds is often deduced 
from 



l&TTcrT^E?{Q^/a) 
L, ' 



1.87 x 10^ 



V 



kms 



days 



(4) 



(sec Gchrz et al. 1980); this follows from equating the to- 
tal power absorbed by the grains in the shell to the total 
re-radiated power. Here L, is the (constant) bolometric lu- 
minosity of the stellar remnant, _R is the grain-nova distance, 
Qc is the Planck mean of the grain cmissivity. Text is the op- 
tical depth due to extinction in the dust shell and 0(roxt) is 
a factor that takes into account extinction within the shell 
(see below). Equation (4) assumes graphitic grains, for which 
Qc oc aTi 

Equation (4) with ^(rext) = 1 follows if (i) there is no 
internal extinction in the dust shell and (ii) the dust shell 
is geometrically thin so that the isothermal grain approxi- 
mation applies. Application of Equation (4) to novae at IR 
maximum leads typically to grain sizes ~ 1 /im, considerably 
larger than even the amax deduced here. For example, in the 
case of V705 Gas, and taking relevant parameters from Mar 
son et al. (1998), we find a(IR max) ~ 1.4/im. 



Infrared spectroscopy of Nova Cassiopeiae 1993 7 



A possible explanation of course is that, as the amax 
deduced here is for t > 253 days, well after IR maximum 
(which occurred around t ~ 100 days in V705 Cas; Mason 
et al. 1998), there was a period of rapid grain destruction fol- 
lowing the grain growth that preceded IR maximum. Indeed, 
application of Equation (4) to V705 Cas for times when 
Text was small suggests that the grain radius initially in- 
creased, and eventually declined. Vaxious scenarios for grain 
destruction have been proposed by Mitchell and cowork- 
ers (Mitchell, Evans & Bode 1983; Mitchell & Evans 1984; 
Mitchell, Evans & Albinson 1986) and it is inevitable that 
grains in nova winds are heavily processed (see Rawlings & 
Evans 2005). 

Alternatively Equation (4) may not be appropriate at 
IR maximum for a dust shell that is optically thick at short 
wavelengths: as we have seen, assumption (ii) above is a good 
approximation in the case of V705 Cas but (i) is not. In this 
case, 

_ 1 - exp[-Text] 

<P(Text) — , 

Text 

which — » 1 as Text — > 0. Applying this correction to V705 Cas 

at IR maximum gives a(IR max) ~ 0.2 ^m for Text — 5.5 (see 
§3.1.3), considerably less than that obtained on the assump- 
tion that the dust shell is optically thin. 

We also note that Equation (4), with 0(Text) = 1, will 
apply whenever L« is constant and the dust shell is optically 
thin at short wavelengths (as is the case after light curve 
recovery), as in V705 Cas at Epochs 1 and 2. Using L, = 
5.6 X 10" L© from §4.2 and V = 850kms"\ Equation (4) 
with 0(Text) = 1 gives 0.5 /um (0.3 /jm) at f = 253 days 
(320 days) , again considerably smaller that that obtained for 
the time of IR maximum. While the situation in V705 Cas is 
complicated by the presence of (at least) two grain types, we 
conclude that the application of Equation (4) to determine 
grain size is inappropriate at infrared maximum in novae 
having optically thick dust shells. 

4.5 The dust composition 

We find that both sets of data arc best fitted by a dust shell 
in which carbon grains are numerically more abundant than 
silicate grains. We do not expect the ACisilicate fraction to 
vary as the dust shell moves away from the stellar remnant, 
and the individual values of the AC and silicate fractions in 
Table 1 are gratifyingly constant; for the purpose of discus- 
sion we average the ACisilicate ratio, and integrate over the 
grain size distribution (see §4.4 below) to get a ~ 0.67 : 0.33 
(± ~ 0.1) carbon:silicatc mix. This translates to a mass 
ratio carbon:silicate = 1.26, significantly lower that the car- 
bon:silicate ratio of ~ 8.5 obtained by Mason et al. (1998). 
However these authors fit a 4-component model to limited 
broadband data, with UIR components at their 'standard' 
wavelengths, and they caution that their derivation of the 
silicate mass may be problematic. We should not therefore 
be surprised by the discrepancy and consider that the result 
obtained here is reasonably robust. 

The co-existence of carbon and silicate dust ('chemical 
dichotomy') in the IR spectra of evolved stars is usually as- 
cribed to the trapping of silicate dust in a disc arising from 
mass-loss during an earlier phase of evolution when the star 
was oxygen-rich, while carbon dust arises from more recent 



or current (carbon-rich) mass-loss (see e.g. Zijlstra et al. 
2001 and references therein). In the case of V705 Cas, UIR 
emission was present in the 10 pm band when silicate was 
weak or absent (Evans et al. 1997). Novae may also have two 
well-separated mass-loss episodes with different chemistries 
to produce their chemical dichotomies, or there may be steep 
abundance gradients in the ejecta, or (as CO formation does 
not go to saturation; see Paper I, Pontefract & Rawlings 
2004, Rawlings & Evans 2005) the nucleation chemistry al- 
lows the simultaneous formation of both carbon and silicate 
grains. In any of these cases, there may be interesting im- 
plications for the thermonuclear runaway, the nature of the 
nova explosion and the pre-dust chemistry. 

While the assumptions leading to Equation (3) may not 
necessarily apply in the case of novae (see Evans & Rawl- 
ings 1994, Rawlings & Evans 2005 for a discussion of this), 
our values for the dust temperature at the inner edge of 
the shell imply condensation temperature Tc — 1060^28 
(1153t4s K) for Epoch 1 (Epoch 2); again the value of Tc is 
not sensitive to and we have assumed l3 — 1. These values 
are of the right order for carbon condensation; however the 
situation in V705 Cas is complicated by the presence of two 
grain types, so the meaning of the 'condensation tempera- 
ture' is not immediately obvious in this case. On the other 
hand, UIR emission and an underlying continuum, but no 
9.7 /xm silicate feature, were present in 1994 May (Paper II), 
implying that the AC condensed before the silicate. 



4.6 The dust features 

4-. 6.1 The silicate features 

In Paper II we suggested that, in addition to the 9.7 /im sil- 
icate feature, there might also be present the corresponding 
18 /xm feature. The IR spectrum, together with the fitted 
model, in the range 7-25 /xm for the two epochs is shown 
in Fig. 3; it seems that the 18 /im feature may indeed be 
present in the IR spectrum of V705 Cas, as may be seen by 
viewing Fig. 3 at grazing incidence. 

As noted in Paper II, Nuth & Hecht (1990) suggested, 
on the basis of laboratory experiments on silicate analogues, 
that the strength of the IS/itm silicate feature relative to 
that of the 9.7 /xm feature increases as the silicate 'ages' 
and is subjected to annealing. Processing of silicates, and its 
effect on the 9.7 /im feature in circumstellar environments, 
has been discussed from an observational point of view by 
Bouwman et al. (2001); however the objects (Herbig Ae/Be 
stars) in this study are young objects, and are not themselves 
dust producers and comparison with the silicates around 
evolved, dust-forming stars, is more relevant from our point 
of view. 

To explore this we also include in Fig. 3 the ISO SWS 
spectrum of the actively dust-forming red supergiant fj, Cep, 
which also displays prominent silicate features (see Tielens 
et al. 1997). We note in particular 

(i) the similarity of the 9.7 /im feature profiles in the two 
objects. The features are narrow, structureless, and peak 
around 9.7 /um, suggestive of amorphous silicate; 

(ii) the relative strength of the 9.7 fw and 18 /ttm features, 
in that the ratio 18/9.7 is much smaller in V705 Cas than 
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Figure 3. The profiles of the siUcate features in V705 Cas for 
the two Epochs; the wavelengths of the 8.1, 8.6, 9.2 and 11.4/im 
UIR features (see Fig. 4b and Table 2) are indicated. The curves 
through the data are the dusty fits; see text for details. Also 
shown for comparison (lower curve) is the ISO SWS spectrum of 
the red supcrgiant fj, Cep. 



in fj, Cep, consistent with the view that the silicate in the 
nova is much 'fresher' than that in Cep. 

While our data lend support to the conclusion of Nuth & 
Hecht (1990), the mid-IR spectra of novae, obtained over 
periods ranging from ~ 20 to 1000 days from outburst 
(and hence over a range of annealing times), suggests that 
the situation is more complex than this. Although V705 Cas 
is now a very weak mid-IR source (the dust was not detected 
with the ISO SWS by Salama et al. 1999), further observa- 
tions in the mid-IR (e.g. with the Spitzer Space Telescope 
(Werner et al. 2004)) will be valuable to monitor both the 
profile of the 9.7 /um feature - which may shift to longer 
wavelengths and display structure as the dust anneals - and 
the development of the relative strengths of the 9.7 um and 
18 /um features. 

4.6.2 The UIR features 

We have subtracted the model fit from the data in the L 
and N bands to reveal the profiles of the UIR emission (see 
Fig. 4a,b). Clearly the resultant UIR profiles in the N band 
depend to some extent on the profile of the warm silicate 
feature (Ossenkopf et al. 1992) but, as we are looking for 
possible changes in the UIR features this does not qualita- 
tively affect the discussion that follows. 

In Fig. 4a the L-band data from Epoch 2 data have been 
normalized to the Epoch 1 data at 3.4 /im and it is evident 
that, while the relative strengths of the 3.28 /^m feature and 
the 3.5 /Ltm 'plateau' are essentially unchanged between the 
two Epochs, the 3.4/3.28 ratio may have changed. 

There is a weak feature at A = 3.522 ± 0.001 /im (see 
Fig 4a). While this is close to the UIR feature at 3.52 /xm 
(see Geballe 1997), the width of the feature in V705 Cas 
(FWHM ~ 0.014 fxm at Epoch 2 and corresponding to an 
expansion velocity ~ 600 km s"'^) suggests that it might be 
a nebular feature rather than a dust feature, which would 
typically be much broader than this (see Table 2 below). 
Further IR spectroscopy of dusty novae in this spectral re- 
gion is desirable to resolve this issue. 



A possible (nebular) identification for this feature is 
Hi 23-6 A3.522 /im, although there is no evidence for Hi 
24-6 A3. 501 pm, which is expected to be of comparable 
strength. An alternative identification, particularly in view 
of the common occurrence of [Fell] lines in the IR spec- 
tra of novae (Rudy et al. 2003; Evans et al. 2003), is [Feii] 
a®Dg/2 — a*F6/2 at A = 3.52367 /im. However there are many 
permitted and forbidden Fell transitions in the wavelength 
range for which we have data and we defer a discussion of 
the emission line spectrum of V705 Cas to a later paper. 

In Fig. 4b the A''-band data from Epoch 2 data have 
been normalized to the Epoch 1 data at 9 /um; the pri- 
mary feature at ~ 8.1 fim has a 'shoulder' at ~ 8.6 fim and 
^ 9.2 /xm. We have fitted the '8.1 /im' feature with 3 gaus- 
sians (see Fig. 5a,b and Table 2). We note from Table 2 that 
the widths of the features are typically ^ 0.3 — 0.6 /im, and 
so these are unlikely to be nebular in origin. Also listed in 
Table 2 are the 'standard' wavelengths and likely identifi- 
cation of 'Class A' UIR features (i.e. those usually seen in 
planetary nebulae and Hll regions) from Geballe (1997). 

There is no evidence that the central wavelengths of 
the two 'shoulder' features (~ 8.72 iim, 9.20 fim) change 
between Epochs 1 and 2; however there seems to be a dis- 
tinct change in the peak wavelength of the primary feature 
between Epochs 1 and 2 (see Fig. 4a and Table 2), from 
8.17 //m to 8.06 /im. This is not an artefact of the subtrac- 
tion of the model from the data, and likely points to a change 
in the character of the carrier. 

As noted in Paper II, the 11.3 /jm feature in V705 Cas 
appears at a longer wavelength than is usual (the slight dif- 
ference between the values in Table 2 and in Paper II is a 
result of the slight diiference in the subtracted continuum). 
We note that there seems to be a significant increase in the 
strength of the 11.3 /xm feature relative to that of the '8.1' 
feature. 



4.6.3 Why are nova UIR features unique? 

It is of interest to consider why nova UIR features are in 
a class of their own (Geballe 1997), and in particular are 
significantly different from those seen in other astrophysical 
environments. 

Laboratory measurements have shown that the peak 
wavelengths of the '3.28', '6.2', '7.7', '8.1' and '11.3' UIR 
features in quenched carbonaceous composite depend on the 
i3Qy'i2Q j-g^^jo^ ^^jjg peak wavelengths increasing linearly with 
the "C/^^C ratio (Wada et al. 2003). However there are 
two reasons why this mechanism can not be operating in 
V705 Cas. First, the "C/^^C ratio in V705 Cas was low, 
^0.2 (Paper I), so isotopic effects will have negligible effect 
on the wavelengths of the UIR features; however this effect 
must clearly be looked for in novae with high ^^C/^'^G ra- 
tio. Second, the peak wavelengths of the nova 3.28 and 3.4 
features are not anomalous, only their relative strengths are 
unusual. If isotopic effects were important the peak wave- 
lengths of all the UIR features would be affected, which was 
not the case in V705 Cas. 

More plausibly, the nature of the UIR carrier in novae is 
governed by the unique environment in which the pre-dust 
and dust chemistry takes place. We consider separately the 
UIR features in the 3 //m and 7-13 /xm windows. 
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Figure 4. The profiles of the UIR features in V705 Cas; solid line, Epoch 1, broken line, Epoch 2. (a) The L-band features; the tick 
mark indicates the feature at 3.522 /um (see text), (b) the AT-band features. Note the shift in wavelength of the 8/um feature. See text for 
details. 



Table 2. Properties of the UIR features in the N-band. 



Epoch 


A 


FWHM 


Flux 


Identification Standard A (//ni) 




(//m) 


(//m) 


(10 Will -) 


lor Class A Tll-is 


1 


8.17 ±0.01 


0.63 ±0.03 


32.1 (±1.4) 


Si— CH3 vibration? 




8.72 ± 0.02 


0.41 ±0.03 


10.8 (±8.1) 


NH2 rock? 




9.20 ±0.02 


0.46 ± 0.04 


6.27 (±0.43) 


? 




11.41 ±0.04 


0.34 ± 0.08 


2.53 (±0.51) 


Aromatic C— H bend 11.3 


2 


8.06 ± 0.06 


0.47 ±0.21 


7.27 (±2.46) 


Si— CH3 vibration? 




8.67 ± 0.08 


0.53 ± 0.20 


4.86 (±1.89) 


NH2 rock? 




9.26 ±0.15 


0.52 ±0.24 


2.34 (±1.29) 


? 




11.39 ±0.12 


0.38 ±0.19 


1.51 (±1.14) 


Aromatic C-H bend 11.3 



4.6.3.1 The 3 /jm window. As is well known, the rel- 
ative strengths of the '3.28' and '3.4' features in novae dif- 
fer substantially from those in other sources such as plane- 
tary nebuleie and post-AGB stars (see e.g. Geballe 1997, van 
Diedenhoven et al. 2004). A possible explanation for this is 
the incorporation of CH2 and CH3 groups in the silicate lat- 
tice, which experimental studies (Grishko & Duley 2002a) 
have shown enhances the 3.4 jim. UIR feature. 

As noted in §4.5, one explanation for the chemical di- 
chotomy in nova dust shells is that CO does not proceed to 
saturation, allowing the simultaneous condensation of car- 
bonaceous and silicate dusts: the chemistry of grain format 
tion and growth is likely to include the formation of carbona- 
ceous grains with silicate inclusions, and vice versa. The 
incorporation of CH2/CH3 groups into the silicate lattice 
in this way is consistent with the simultaneous presence of 
C-bearing and Si/O-bearing dust, and underlines the irrel- 
evance of the C;0 ratio to the nature of the condensate in 
nova winds. 

Furthermore, (see §4.4), the grain size distribution indi- 
cates that larger grains may have swept up the smaller. Un- 
less there were (for example) substantial species-dependent 
grain charge effects it is unlikely that larger grains would 
discriminate between one grain type and another in the 
sweeping-up process; this would also give rise to the kind 
of composite grain hinted at by the 3 /xm UIR features. 



4.6.3.2 The 7-13 /iim window. Additional support for 
the hypothesis that CH2 /CH3 inclusions in silicates could be 
significant in novae comes from the '8.1' feature in V705 Cas. 
A feature at this wavelength is also seen in laboratory spec- 
tra, and may be identified with Si— CH3 vibration (W. W. 
Duley, private communication) . Indeed Duley (private com- 
munication) has suggested that the shift in the peak wave- 
length of the 'S.r feature in V705 Cas between Epochs 1 
and 2 (see Fig. 4 and Table 2) may be the result of a change 
in the co-ordination of Si— CH3 groups. 

Furthermore, it is likely that the UIR carrier in 
V705 Cas (and indeed in other novae) is heavily contam- 
inated by other species, notably O or N, particularly the 
latter given its extreme overabundance (~ 100; see sum- 
mary in Gchrz ot al., 1998) relative to solar in nova winds. 
Grishko & Duley (2002b) have investigated experimentally 
the IR spectra of HAC prepared in the presence of various 
contaminants, and find that HAC deposited in the presence 
of N-bcaring molecules (particularly NH3) displays several 
features not present in 'normal' UIR carriers. In particu- 
lar, there is a feature at 8.56 /im (which Grishko & Duley 
attribute to NH2 rock), and a broad (unidentified) feature 
centred at ~ 9.1 /um with width 0.4 yum, which are close 
to the features on the 'shoulder' of the 8.1 /im feature in 
V705 Cas (see Fig. 5a,b and Table 2). We tentatively sug- 
gest that the carrier of the UIR features in V705 Cas had a 
significant nitrogen component. 
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5 CONCLUSIONS 

Wc have used the dusty code to reassess the properties 
of the optically thick dust shell around the classical nova 
V705 Cas. We find 

(i) the dust shell is a ~ 2 : 1 (by number) mix of amor- 
phous carbon and 'warm silicate'; 

(ii) the grain size distribution n{a) da oc da was flat, 
with q ~ 2.3; this may indicate that the larger grains swept 
up the smaller grains; 

(iii) the maximum grain size in the distribution was 
Omax — 0.06 nm, considerably less than that deduced from 
previous studies (Shore et al. 1994, Paper II, Gehrz et al. 
1998); 

(iv) the dust-bearing shell was, like the CO-bearing shell, 
geometrically thin, such that AR/R ~ 0.1; 

(v) the condensation temperature of the dust was ~ 
1100 K, although in V705 Cas the situation is complicated 
by the presence of two condensates; 

(vi) structure in the 8-9 fiin UIR complex may be due to 
nitrogenation of the UIR carrier; 

(vii) the strength of the 3.4 /mi feature relative to the 
3.28 /Ltm feature, commonly seen in novae, as well as the 
8.1 /um feature, may be due to the incorporation of CHg and 
CHa groups in the silicate matrix; 

(viii) the silicate dust at least weis freshly condensed, 
judging by the weakness of the 18 /um silicate feature; 

(ix) the extrapolated angular diameter of the dust shell 
is such that the resolved remnant seen by Diaz et al. (2001) 
was very likely due to line emission rather than dust. 
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